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.  .'TRACT 

Geostrophic  relative  vorticity  is  comruted  along  three  latitude 
circles  generally  encompassing  the  United  States,  using  a  Fourier  analy- 
sis of  the  500-mb  height  field  for  the  period  3-6  November  1961. 

A  statistical  comparison  of  vorticity  obtained  in  this  manner  with 
it*s  square  grid  computation  is  made. 

The  writers  are  deeply  indebted  to  Professor  F.  L.  Martin  of  the 
United  States  Naval  Postgraduate  School  for  his  suggestion  of  the  topic 
and  his  continued  help  throughout  the  investigation  and  during  the 
preparation  of  this  paper. 
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1  .   Introduction. 

The  advent  of  the  electronic  computer  coupled  with  increased  data 
coverage  has  removed  the  primary  obstacles  to  the  use  of  harmonic  analy- 
sis as  a  useful  instrument  in  atmospheric  research. 

Use  is  made  of  a  Fourier  analysis  of  the  500-mb  height  field.  A 
computer  program  devised  by  G,   Arnason  and  M.  Reese. and  made  available 
by  Fleet  Numerical  Weather  Facility'  was  used  to  obtain  the  harmonic 
analysis  around  latitude  circles  at  distances  of  1  through  27  grid  lenrths 
d  =  3#1  km  from  the  pole.  An  equation  for  relative  vorticity  using  the 
results  of  the  spherical  harmonic  analysis  is  derived  and  tested  statisti- 
cally against  the  more  usual  numerical  computation  of  <~eostropbic  relative 
vorticity  made  using  a  square  grid. 

The  problem  undertaken  was  to  compute  relative  vorticity  at  500  mb 
using  the  geostrophic  vorticity  equation  in  its  spherical  coordinate  form 
[1,  P.  358], 

l~^(~\   +-  ^    )       ^  Tan   *  H  (1) 

In  equation  (1  ),  z  was  replaced  by  its  expression  in  terms  of  a  Fourier 
analysis 

F.~    ?f   [[cnii^(fiA-en)]  ~    f   f  /*  (2) 

Here  z  is  the  zonally  averaged  contour  height,  whereas  z*'  is  the  perturba- 
tion quantity  which  remains  when  the  z  field  is  subtracted  from  the  actual 
height  field.  The  computer  printout  included  maps  of  each  of  these  three 
fields  at  500  mb  for  the  period  under  investigation,  3  through  6  November 
1961.  Geostrouhic  relative  vorticity  was  computed  by  the  above  mentioned 
technique  for  17  longitudes  along  each  of  three  latitude  circles,  generally 

1 
Henceforth  referred  to  as  FNWF 


encompassing  the  north-south  extent  of  the  United  States.  This  gave  a 
total  sample  size  of  204  computed  relative  vorticities  which  were  then 
correlated  with  square-grid  computations  of  relative  vorticities  obtained 
from  FT-IVJF  charts  of  absolute  vorticity  for  exactly  the  same  set  of  points 
and  times. 

2.  The  geostrophic  relative  vorticity  in  terms  of  a  harmonically  analyzed 
contour  field. 

The  adequacy  of  the  Fourier  representation  of  meteorological  para- 
meters will  not  be  discussed  here;  the  interested  reader  is  referred  to 

L2  j . 

It  was  previously  mentioned  that  use  was  made  of  the  FNUF  computer 
program  for  a.  Fourier  analysis  of  the  500  mb  contour  field.  A  brief  de- 
scription of  the  theory  and  numerical  evaluation  of  the  Fourier  coeffi- 
cients is  given  in  [3]   and  only  the  applicable  portions  of  the  program 
will  be  mentioned  here. 

The  zonal  mean  height  'z  is  defined  as 

2  TT 

?  =.^f  I    O  (3) 

and  its  zonal  perturbation,  in  accordance  with  equation  (2),  is 

l^-Z-l  CO 

Expanding  z  in  a  Fourier  series  one  obtains  at  a  given  latitude 

Z  ~  ?  1-  LEA,   '  h*  +~  ^n  ec'5  V^J  ^5) 

where  A  and  B  are  the  coefficients  of  the  nth  harmonic  functions  in  the  • 

r,        n  — • 

Fourier  expansion  of  z,  and  M  is  the  total  number  of  harmonic  waves  neces- 
sary to  resolve  the  contour  field.  Replacing  A^  and  B^  in  (5)  by  an  am- 
plitude factor  C  and  a  phase  angle  ©  .  (5)  mav  be  rewritten  as 

1  \  n 


:.     +-  ^  Ln  r/-v  . //;^  -  :  (6) 


w>>ere 

C,^  lA^-8;  and   tan   ©„■=  -  Bn/*„ 

Equation  (6)  has  been  programmed  for  the  CDC-16G4.  digital  computer  by 
FNWF,  and  the  program  includes:  a  contour  grid  map  of  the  zonal  pertur- 
bation field  Z'f,  and  a  tabular  list  of  the  Fourier  parameters  Cn,  On  and 
If  for  each  of  27  latitude  circles .  The  last  three  parameters  are  actually- 
computed  in  relation  to  a  zonal  band  of  width  equal  to  the  grid  mesh  d  of 
the  FMJF  map  grid,  but  may  be  considered  applicable  at  the  center  of  the 
band . 

Using  spherical  coordinates  A  and  (b  (longitude  and  latitude,  respec- 
tively) the  derivative  operators 

i-  =  J~  i         and'      i-:   ■'  ■  >  - 

apply,  and  we  have 

2,  U,   .5    .A         (7) 


:   ,   /  -L  •  15.  f  i>.  "j  f  f 


since 


If  now  z  is  replaced  in  (7)  by  its  expansion  in  a  Fourier  series  as  given 
in  equation  (2)  we  obtain 

+-  |-r  [_}   (1  5"<  w^-  fc,Q(  f  f  Uj  cuic  ;  f  (8) 

Remembering  that  C  and  :]    are  independent  of  longitude  })   but  not  of 


latitude  .[  ,  we  may  perform  the  indicated  differentiations  and  combine 
terms,  with  the  result  that  J*  in  spherical  coordinates  becomes 
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-[l^s(„*         :,  *,  ^tC|   1;  j   ■    ',«:.., 


Equation  (9 )  may  be  shown  to  have  the  form 


(9) 
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V  I 
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-|^)"jJCn  im(^-\) 


(10) 
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equation  (K  )  may  be  approximated  by 


-  »'i 
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L  ^ 
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(11) 


4   M 

Here  the  superscript  bar  denotes  an  average  with  respect  to  n  over  all 
significant  vraves  n  =  1,  ..„..,  M. 

The  FNWF  maps  are  stereographic  maps  true  at  latitude  60N,  and  else- 
where the  grid  mesh,  d  =  381  km,  is  related  to  the  true  distance  dT  by  the 

relationship  d  =  md_ ,  where  m  -  (1  +  sin  60)/(l  +  sin  h)  and  is  the  so- 

i  r 

called  map  factor. 

The  partial  derivatives  in  (11  )  may  be  evaluated,  by  centered  finite 
differences.  The  numerical  values  of  C  n  and  <^,x  were  simply  taken  from 


the  tabular  list  of  computer  output  of  the  harmonic  analysis  at  the  vari- 
ous latitude  circles.  For  example,  )  '  f      may  be  written  using 
finite  differences  between  adjacent  latitude  circles  (for  which  QA$      d/m  ) 
in  the  form 

The  latitude  circles  of  the  Fourier  analysis  employed  in  the  compita- 
tion  of  1^  corresponded  to  N  =  12,  15,  and  18  grid  mesh  distances  from  the 
pole,  which  in  turn  corresponded  to  the  latitude  circles:  A.8.9II,  3S.6N, 
and  28, 3N  (see  fig,  1).   However,  in  the  computation  of  first  and  second 
derivatives  of  the  nth  Fourier  wave  parameters,  these  parameters  at  lati- 
tude circles  N-1  and  N+1  are  also  needed. 

The  term  at/   f\     is  more  tractible  in  the  present  case  if  the  finite- 
difference  substitution  r X  }\    -  La/m,    is  made;  thus,  the  centered  deriva- 
tive °*l/&h     in  finite  differences,  assumes  the  form 


'tiC&-~<)  '     (12) 


where  "f  -  $  ■■:- ■    j   is  the  radius  of  the  latitude  circle,  and  the  subscript 
i-i-1  is  one  grid  mesh  eastward  of  the  ith  uoint  along  the  latitude  circle. 

Including  the  form  of  s ,  V  . -'  riven  in  (12),  equation  (11  )  may  now 
be  written  in  the  more  compact  form 

A  -  V  "  i  ^c:  t    •  ^.2?*-  »'..  (.*£,-  **->)         (13) 


where 
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The  averaging  with  respect  to  n  in  the  expressions  for  W.  and  W?  is 
assumed  to  be  valid,  according  to  the  mean  value  theorem.  Actually  how- 
ever the  average  was  approximated  by  a  simple  arithmetic  mean  of  the  con- 
tributions for  n  =  1,  . ...,  M,  Thus  it  was  necessary  to  compute  contri- 
butions for  each  n  and  divide  the  result  by  the  number  of  waves* 

Van  Mieghem  j^2J  points  out  that  waves  of  space  number  higher  than  12 
to  17  contribute  little  to  the  variance  of  the  contour-height  field.  Con- 
sequently, the  Forrrier  analysis  was  truncated  at  M  =  12  at  the  two  higher 
latitude  circles,  although  wave  numbers  1  through  15  were  employed  at  the 
lowest  latitude. 

The  coefficients  W.  and  Wz  are  constants  for  a  particular  time  and 
latitude  circle,  and  can  be  computed  directly  from  the  tabulated  harmonic 
data.  These  values  are  listed  in  table  1 ■   However,  in  applying  equation 
(13)  it  was  necessary  to  interpolate  values  of  z";;  from  the  perturbation 
map  at  the  points  indicated  in  fig.  1  along  the  specified  latitude  circles, 
T^ese  points  were  arbitrarily  chosen  at  integral  multiples  of  five  degrees 
of  longitude,  encompassing  the  west-east  extent  of  the  United  States, 
Therefore,  for  ease  in  selecting  data  points,  a  grid  ovcrlajr  was  con- 
structed presenting  all  latitude  circles  Nf=11  ,*..*,  19,  drawn  so  as  to  in- 
corporate the  map  factor. 


•1 
At  the  two  higher  latitudes,  W,  was  negative  ensuring  that  lo\/s 

corresponded  to  positive  relative  vorticity.  At  the  lowest  latitude 
(N=18),  it  was  necessary  to  use  M=15  in  order  to  ensure  that  this  condi- 
tion was  met., 


:   IS  1:  Values  of  '/,  and  T.J2  to  be  applied  in  equation  (;  ).   All  values 

—i  _  j        1 

are  in  units  of  1 0   sec  "  f  t  . 

r    ]   -  -— ^_.„_-. p— 

,_ , ! r , 1 

DATE  N  =  12         N  =  15         N  =  18 


3  Nov.  1961 


W  =  -.U67     W  =  -,1635     W  =  -.2696 


W  =  -.0079   !  W  =  -.0370     W  =  -.0303 

j 


W 

= 

-.1490 

w 

= 

-.0617 

w 

= 

-80420 

k 

Nov. 

1961 

w 

= 

-.0311 

w 

= 

.0316 

w 

= 

.0138 

5  Nov.  1961 


6  Nov%  1961 


W  =  -.1355  W  =  -.1013  W  =  -.1819 

W  =  -.0398  W  =  -.0230  W  =  .1682 

T"  =  -.0551  "=~-7u60  w  =  -.3510  \ 

W  =  .0757  W  =  .0349  W  =  -.7610 

Note  that  the  term  (2  '\,  s>  Z$)/d      in  both  (11  )  and  (13)  also  varies 
from  point  to  noint  around  any  latitude  circle*  However,  the  variation  of 
this  term  was  so  small  that  the  mean  zonal  wind  for  the  latitude  was  used* 
The  mean  zonal  x;ind  TJ*  was  computed  using  the  latitudinal  gradient  of 
zonal  mean  height  z,  that  is  from 

again  using  finite  differences , 

Along  each  of  the  three  latitude  circles,  17  values  of  W,  z*  and 
W  (z*  -z*,  )  were  obtained  for  each  of  the  four  days,  living  a  total  of 
204.  separate  values  of  these  quantities. 
3.  Statistical  results. 

For  precisely  the  same  set  of  points  and  days  as  described  above, 
corresponding  values  of  J  3  were  obtained  from  grid  maps  of  geostrophic 
relative  vorticity  made  available  by  FNIJF.  Their  maps  essentially  pre- 
sented the  field  of  j«  computed  using  their  square-grid  mesh  according 
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to  the  usual  formula 


U  =  —v1-  /z  -0.  <1i) 


where  z  is  the  four-point  space-mean  contour  height  and  d  =  381  km  at 
latitude  60N, 

Since  the  averaging  process  in  equation  (11  )  is  somewhat  approximate, 
it  was  felt  that  better  verification  would  result  if  the  relationship  be- 
tween J  wf=  and  W  zf   and  W„  (a£  -zf.  )  v/as  written 

Here  Ao  ,  A, ,  and  A?  are  regression  coefficients'  to  be  determined  by  a 
least  squares  technique  applied  to  the  204,  simultaneous  values  of  W,  zf 
and  W,  (a*  -a*  )«  Because  of  occasional  large  gradients  of  3_  ,„  over  5° 
longitude  intervals,  the  values  of  /  ^  were  smoothed  by  the  three-point 
binomial  running  mean  formula 

where  X;_  is  any  value  of  j  a    interpolated  from  the  FNl-IF  vorti city-grid 
map  at  each  point.  The  regression  equation  actually  involves  the  correla- 
tion of  the  smoothed  )  values  with  the  independent  variables  already 
noted* 

Using  the  BIMD  06  statistical  program  as  adapted  to  the  CDC-1604. 
computer,  it  was  necessary  to  type  204.  IBM  cards  containing  the  smoothed 
J-./VkVp  values  and  the  simultaneous  values  of  the  independent  variables  of 
the  right  side  of  equation  (15)-  Best-fitting  values  of  the  regression 
coefficients  were 

•1 

'These  regression  coefficients  are  non-dimensional  since  the  units 
of  W,  z?   and  W„(a*  -z:"",)  are  in  sec"^  . 
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A0  =  -8.1719,  &|  =  O.4535,  and  A2  =  0,0624- 

The  standard  error  of  estimate'  and  the  multiple  linear  correlation  co- 
efficient for  the  regression  equation  were  20.3523  and  0.7632  respectively. 

Panofsky  and  Brier  [4.,  p.  11 3J  present  an  F  test  for  the  analysis  of 
the  variance  explained  by  the  regression  equation.  The  F  value  may  be  de- 
fined by  the  ratio  of  the  explained  variance  to  the  residual  variance 


relative  to  the  recession  equation.  This  formula  is 


f-    *.yxt\v    .  (16) 


With  p  =  2,  the  number  of  independent  variables,  and^]  =  204,  the  number 
of  random  samples,  the  value  of  F  from  equation  (16)  turns  out  to  be 
139.33.  however,  assuming  normally  distributed  independent  and  dependent 
variables,  the  critical  F  at  the  0.99  significance  level  with  2  and  201 
degrees  of  freedom  is  /+.71.  Thus  the  regression  equation  (15)  is  signi- 
ficant at  a  very  high  confidence  level. 

The  partial  correlation  coefficients  of  W  zf  and  W„  (z*  -z*  )  were 
0.7620  and  0.0382  respectively,  indicating  very  little  contribution  to 
the  regression  equation  for  the  latter  variable .  An  intuitive  explana- 
tion for  the  insignificant  contribution  of  W, (z*  -z*.  )  mav  be  obtained 
from  fig,  2,  The  differences,  z?\  -z-\  ,  taken  over  the  pressure  centers 
would  be  zero,  whereas  in  actuali'y  the  largest  absolute  values  of  j 
would  be  found  at  these  points,   likewise  for  a  sinusoidal  contour  wave 
pattern  maximum  magnitudes  of  z*   -z::"  would  exist  near  contour  inflection 
points,  where  d  1 

The  possibility  of  "automatic"  correlation,  viz,  correlation  between 
the  two  independent  variables,  was  examined  and  none  was  found  to  exist. 

All  values  of  £«.  are  in  units  x  10  °  sec""1  . 
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Fir.  2.  Distribution  of  j^  along  a  latitude  circle  between  contour  maxi- 
mum  and  minimum,  in  the  case  of  a  sinusoidal,  contour  wave. 

The  individual  latitude  circles  were  statistically  examined  in  a  like 
manner.  The  two  northernmost  latitude  circles  gave  a  multiple  linear  cor- 
relation coefficient  of  O.84.,  whereas  there  is  no  significant  correlation 
between  JPA/M/,  and  the  two  independent  variables  W,  zf  and  W,  (z;;"  -zx'  )  at 
latitude  2S.3N.  This  inconsistency  at  low  latitudes  is  attributable  to 
"noise"  in  the  Fourier  representation  of  the  contour  field.  The  "noise" 
is  due  to  the  small  variance  of  the  height  field  at  this  latitude,  coupled 
with  the  relatively  large  RMS  error  possible  in  the  500  mb  analysis  in  the 
subtropics  during  the  period  of  study.  The  use  o^  second  derivatives  of 
the  Fourier  parameters  also  increases  the  amount  of  "noise"  generated  in 
the  determination  of  J  ,  , 
4-.  Conclusions, 

Theoretically  there  should  have  been  a  large  positive  multiple  linear 
correlation  between  J     and  W  z5,-  and  W„  (z*  -z.  ,)»  This  Droved  to  be  the 
case  at  latitudes  4-2. 9N  and  38. 6N  where  the  contour  field  was  well  resolved 
by  twelve  harmonic  waves.  Fowever,  at  the  lowest  latitude  circle,  28. 3N, 
the  several  sources  of  errors  listed  at  the  end  of  section  3  could,  and 
probably  did,  influence  the  low  correlations  that  were  obtained.  These 
were  basically  due  to  the  lack  of  resolution  of  the  contour  field  by  the 
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Fourier  analysis.  However,  the  sane  reasons  would  cause  the  square-grid 
computation  of  \ -,    to  be  subject  to  large  errors. 

Based  upon  the  results  at  the  two  higher  latitudes,  it  seems  apparent 
that  the  use  of  a  Fourier  analysis  of  a  meteorological  field  presents  a 
very  valuable  tool  for  atmospheric  research,,  With  the  advent  of  the 
electronic  computer  and  machine -analyzed  contour  charts,  the  use  of  har- 
monic analysis  is  an  effective  and  economically  feasible  approach  to  the 
computation  of  synoptic  aids  for  the  practising  meteorologist. 
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